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I. INTRODUCTION 

Pion production in nucleon-nucleon collisions has been 
given increased attention due to the recent availability 
of precise near-threshold data for pp —> ppir ', pp — > dn + 
and pp — ► pnir + . These reactions involve a large momen- 
tum transfer, but near threshold only a few partial waves 
can contribute. Here the study of the short-range fea- 
tures of the NN interaction is greatly facilitated because 
it is possible to experimentally separate the contributions 
from different angular momentum states. In recent stud- 
ies it was found that within 10-20 MeV of threshold, 
i.e. below 300 MeV, only Ss for pp —* ppn°, [QH or 
at most two partial waves (Ss and Sp) for pp — ► dir + and 
pp — > pnn + P-[icfl were important, and direct compar- 
isons with simple theoretical approaches can be made. 
The expected dominance of Ss contributions (Inn = 0, 
Ltt = 0) for the final states was seen, but at the same 
time it was found that for pp — > ppir traditional models 
under- predicted the cross sections by factors of 3-5. 

At higher energies the reaction pp — > pnn + tends to be 



dominated by the delta resonance. The present 

study in the near threshold region will provide informa- 
tion on polarization observables in a region still subject 
to relatively parameter-free microscopic calculations. Its 
major motivation is the hope and expectation to reduce 
the theoretical uncertainties for a process fundamental 
in nuclear physics. Additional references to prior exper- 
imental and theoretical work can be found in fL(i|Js|]. 

Much effort has gone into attempts to explain 
pp^ppir , probably the theoretically most difficult of 
the three NN — ► NN-7T branches. Modern chiral perturba- 
tion work on this reaction has been discussed recently in 
jl4],[l5| , but this approach has not yet explained the data. 
It has become apparent that for pp — > pnir + the first or- 
der terms dominate and delta and complicating heavy 
meson exchange contribution are relatively small making 
this reaction theoretically important. At this time there 
is no refereed publication on chiral perturbation calcula- 
tion; however, very recently a preprint has become avail- 
able that reports xpT calculations for pp — > pnn + and 
pp — > dn + . These xpT calculations under-predict near 
threshold data by factors of 2 to 3, but with less room 
to adjust poorly known parameters than for pp — ► ppir 
Insf . The varying theoretical approaches proposed for 
the inclusion of off-shell contributions, "heavy meson" 
exchange, and contributions involving the A isobar are 
still under discussion |l7j . Divergent results suggest that 
conclusions about the adequacy of a theoretical pion pro- 
duction model cannot be based on agreement for just 
a single pion production branch. A satisfactory model 
must succeed in 

explaining all NN — > NN7r branches. 

Promising attempts to simultaneously describe data 
from more than one reaction have been made using the 



Hamiltonian approach . A relativistic phenomenologi- 
cal model was offered by Shyam and Mosel fjs|] . Their 
approach naturally included higher angular momenta 
and the exchange of heavier mesons and was able to fit 
pp —>■ pnir + and pp — > ppir cross sections over a wide 
energy range. This model should also predict analyzing 
powers; however, such predictions have not been pub- 
lished to date. More complete microscopic calculations 
in the meson exchange picture, involving angular mo- 
mentum states as large as L T = 2 are now being pub- 
lished by the Julich group p9[-pl|. Their approach of 
p0| , pH is largely free of adjustable parameters. It gives a 
very good reproduction of the published pp — > ppir and 
pp — > pnir + data. However, it has only qualitative suc- 
cess in reproducing the new spin correlations observables 
p2f for pp — > ppTT . The spin correlation data reported 
here for pp — > pnit + present a different test of the pre- 
dictive power of this approach. 

Precision measurements of pion production near 
threshold became possible with the advent of cooled 
beams and the technology developed at IUCF and other 
facilities around the world to use a stored cooled beam 
on a windowless internal target. In the present experi- 
ment we used a polarized beam on a polarized internal 
target to measure the pp — > pnir + reaction at 325, 350, 
375 and 400 MeV. Bombarding energies for the present 
experiment were chosen to cover the range over which 
some higher partial waves gradually become significant. 
This is the first investigation of spin correlation coeffi- 
cients in this reaction near threshold. We have measured 
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kinematic variables and the observed range of the of the 
polar angles. We also determined the differential analyz- 
ing powers A y (0). 



II. EXPERIMENT 

The present experiment was conducted at the Cooler 
ring of the Indiana University Cyclotron Facility (IUCF). 
A beam of polarized protons was accelerated to 197 MeV 
and stack-injected into the ring. There it was accumu- 
lated for 2-3 minutes, yielding orbiting beam currents of 
100-300 /iA. The stored beam was then accelerated in 
the ring to full energy. After about 5-10 minutes of data 
taking the remaining beam was discarded and the cycle 
was repeated with reversed beam spin. The beam polar- 
ization was vertical, alternately up or down. 

The target consisted of 74 ± 4% polarized atomic hy- 
drogen in a storage cell, where it was maintained by an 
atomic beam p3[ . The open-ended cell, a movable cylin- 
der of 12 mm diameter with a 25 (im thick aluminum wall, 
minimized the background induced by the beam halo in 
the cell wall. Depolarization due to wall collisions was 
minimized by coating the cell with Teflon. The target gas 
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density distribution along the storage cell was of trian- 
gular shape with a spread of ±12.5 cm about the center. 
The thickness of the target was « 10 13 atoms/cm 2 . The 
target polarization was flipped every 2 sec, pointing in 
sequence up, down, left, right. Each change of direction 
took less than 10 ms. 

Elastic proton-proton scattering events were observed 
concurrently with pion production in order to deter- 
mine the luminosity and beam and target polarization. 
Four plastic scintillators at azimuthal angles of ±45° and 
±135° detected coincident protons elastically scattered 
near 8\ a b — 45°. The product PQ of the beam (P) and 
target (Q) polarization was determined from the spin cor- 
relation A xx — A yy , which is large and well known for pp 
scattering. Typically, we saw PQ « 0.5, but for the 350 
MeV run PQ dropped to 0.34. The time-integrated lu- 
minosity was deduced from the known pp elastic cross 
section. 




FIG. 1. Apparatus for the experiment: WC1 and WC2 are 
wire chambers. E and K are segmented plastic scintillators 
stacks that determine the energy of the charged reaction prod- 
ucts. V is the veto detector. F is a timing detector. X is one 
of the four detectors for the elastic pp scattering monitor. 

The detector stack is shown in Fig. 1. The trajecto- 
ries of the reaction pions and protons were determined 
by wire chambers WC1 and WC2, spaced 23 cm apart, 
each consisting of two wire planes. The angular accuracy 
of the deduced trajectories was limited by multiple scat- 
tering of the ejectiles and by the resolution of the wire 
chambers. The major contribution to the multiple scat- 
tering came from the 1.5 mm thick plastic scintillator de- 
tector (F) which served as timing detector. We estimate 
the angular uncertainty for the trajectories for protons 
and pions of about a = 0.5° and a = 1.0°, respectively. 
The energy of the charged particles was measured by a 
stack of segmented plastic scintillators of total thickness 
of 25.4 cm. The energy resolution of the stack was about 
a = 3.5%. Events are vetoed if a charged particle reaches 
the detector V behind the K-detector. 



Various types of events were analyzed for the present 
experiment. The ]3p — > pn7r + candidates had to have si- 
multaneous counts in one or more F segments and two 
or more E segments. They were rejected if there was a 
count in the veto detector V. The veto detector effec- 
tively eliminated (or else selected) very energetic (elas- 
tic) projectiles. For the 400 MeV run it also rejected the 
most energetic reaction protons and pions, complicating 
the interpretation of this run. In addition to the pion 
trigger there were a "monitor trigger," to study pp elas- 
tic scattering, and a "single prong" trigger for diagnostic 
purposes. 

More details of the experimental set-up are given in 




III. ANALYSIS 

Coincident protons and pions were identified by their 
deposited energy and their time of flight between the E 
and F detectors. Software cuts were used to further elim- 
inate background. We employed a particle identification 
cut, a cut based on the trajectory origin, and a missing- 
mass cut. Fig. 2 shows the distribution of the missing 
mass reconstructed from the measured momenta of the 
proton and pion. 
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FIG. 2. Distributions of the calculated mass m x of the un- 
detected particle at 325 MeV bombarding energy, shown for 
the four combinations of vertical beam and target polariza- 
tion. A sharp peak (» 3.5 MeV/c 2 FWHM) is seen at the 
neutron rest mass. The shaded region indicates the assumed 
background distribution. 
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The events of interest are those near the neutron 
mass. The software cuts and the coincidence require- 
ments greatly reduce background, but some of it remains 
under the neutron mass peak. This amounts to 3 to 8% of 
the peak area and primarily comes from pions produced 
in the walls of the target cell. To assess the background 
shape we took data with pure N2 gas in the target cell. 
We deduced the background under the missing mass peak 
by using a background shape derived from the N2 miss- 
ing mass spectrum. This empirical shape was normalized 
to the visible background continuum outside the neutron 
peak. To estimate the error from background corrections 
we varied the background subtraction by ±25% . The 
effect on the final results was found to be smaller than 
the statistical errors. At 325 MeV the resolution of the 
neutron missing mass peak was a = 1.6 MeV. 

A Monte Carlo simulation of the experiment was used 
to determine various limiting effects of the apparatus and 
to derive corresponding corrections. The code contained 
the detailed geometry of the detector systems and the 
target density distribution. It included the effects of mul- 
tiple scattering of the charged particles, pion decay, loss 
of energy by the charged particles, and the energy res- 
olutions of the scintillators. The simulation provided a 
guideline for the expected angular and energy distribu- 
tions of pions and protons. This was important at 400 
MeV, because about 0.2% of the coincident protons and 
14.4% of the pions, intercepted by the detectors, did not 
stop in the scintillators. In the final analysis we used only 
those events where the proton and the pion energies were 
measured accurately enough to produce a missing mass 
value within ±3.5 MeV of the neutron mass. The simu- 
lation was also needed to study the acceptance of the de- 
tector system, which did not cover the entire phase space 
of the exit channel. Very forward events were suppressed 
by the central hole in the detector which was needed to 
accommodate the beam pipe for the circulating beam. 
Many large-angle events are lost because of the finite de- 
tector size. Coincident events are detected only if the 
polar angle of the proton is larger than 6° and if the pion 
angle is 6° < 6\ a b < 32°. The event acceptance for p7r + 
coincidences was 21.5% at 325 MeV, 19.4% at 350 MeV, 
18.2% at 375 MeV, and 15.0% at 400 MeV. The resulting 
loss in count rate is significant; however, since target and 
projectile are identical particles we sample a phase space 
twice as large as the acceptance suggests. 

In the cm. system low energy pions are seen for all an- 
gles, but energetic pions are not sampled for polar angles 
near 6 c . m . = 90°. Near threshold the analyzing power for 
pions in the pn7r + final state varies about linearly with 
pion momentum, independent of bombarding energy (jjj. 
Therefore A y data for cm. angles between 65° and 115° 
not only have low statistics, but they also require (model 
dependent) corrections, which results in large error bars. 
The integrated spin dependent cross section ratios S, D, 
and A y reported here lack data near 6> c . m . = 90°, and 



therefore only approximately represent the full integral 
over # c .m.- Nevertheless, accurate comparisons with the- 
ory are possible. They require a corresponding integra- 
tion range for the theoretical calculation, i.e., the restric- 
tion of pion lab angles to the range 6° to 32°. 

IV. RESULTS AND DISCUSSION 

Fig. 2 compares four missing mass spectra with the 
same integrated luminosity for the four spin-spin ori- 
entations with vertical beam and target polarization. 
The spin-dependent cross section Aa' T (defined as the 
observed cross sections with the spins of the colliding 
protons opposite minus the observed cross section with 
spins parallel) can be calculated by using the cross ratio 
method 

Aa' T 2 1 - VR , 

^"unpol 

PQ 1 + VR 

where R = (Y uu Ydd/Y u dYd u ). Y mn denotes the yields 
for the four spin combinations and a unp oi is the unpo- 
larized cross section for the measured range of the polar 
angle 0. Because of the limited range of the detector 
the integrated cross section ratio ACT^/cr unpo i will differ 
from the standard ratio Aerx/ctotai if the spin correla- 
tion coefficients have a strong 9 dependence. The recent 
Hamiltonian calculation by Hanhart et al. |^l| suggests 
relatively small changes for our geometry (see Fig. 4). 

For the present experiment spin-dependent differential 
cross section as a function of the pion azimuthal angle 
fa can be expressed as pq] 

a(fa) = cr tot [l + (p y Ay - q y Aj) cos fa 

/ -f" Ayy A XX Ayy . 

+Py1y( Y~ — cos2fa) 

+terms involving p x ,Pz,<lx,<lz}, (2) 

where Ay and Aj are respectively the beam and target 
analyzing powers, and p X) Pz)qxjqz are very small. The 
values A nn are the spin correlation coefficients and p;(qi) 
are the beam (target) polarization components. In this 
equation the observables Ay, Ay and A nn are integrals 
over all kinematic variables of the proton, the neutron 
and the pion except over the pion azimuthal angle 
Spin correlation observables are deduced by sorting the 
data as a function fa for a specific combination of beam 
and target spin. We show the ^-dependence of the mea- 
sured ratios in Fig. 3. 

Except for the scales, the presentation is identical to 
that in ref. p2| for pp — > pp7r°. The columns represent 
the four bombarding energies, and the rows the differ- 
ent combinations Rj for the spin-dependent yields. The 
definitions of the Rj are given in Fig. 3. To obtain the 
plotted quantities, we calculate Wi = (Ri — l)/(Ri + 1) 
and divide by ^PQ for the first three rows and by PQ 
for the last two two rows. 
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FIG. 3. Asymmetries for different spin combinations Ri (listed on the left) as a function of the azimuthal angle of the 
7r + . The solid curves represent a least-square fit using the theoretical <f>-n dependence (listed on the right), varying &b, &t, 
S ee (A xx + Ayy), and D = (A xx - A yy ). 
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The first three rows are directly related to the beam 
and target analyzing powers, whereas the last two rows 
reflect the spin correlation coefficient combinations S = 
(A xx + A yy ), and D = (A xx — A yy ). In this figure as and 
&t are related to the beam and target analyzing powers 
A^ and Aj by a B = A^P/VPQ, and a T = AjQ/v^Q. 
The solid lines in Fig. 3 represent a least-square fit using 
the (j>Tr dependence from Eq.(2) (listed to the right of 
Fig. 3), and varying the four parameters as, aT, S, and 
D. 



Aa T /a„, = -(A„+A yy ) = -S 
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those in Fig. 3, are corrected for background contribu- 
tions. The values for the integrated pp — > pnir + spin 
correlation coefficients in Fig. 4 are in remarkably good 
agreement with recent Hamiltonian predictions following 
ref. |§. 

In our experiment beam and target analyzing powers 
are related by Ay (6) = —AJ(tt—6) and, therefore, permit 
a study of the back angles for A y (0 c . m .) with the present 
apparatus. Fig. 5 shows the observed 9 c . m . dependence 
of the analyzing powers A y (9 c . m .). As mentioned above, 
the limited size of the detector reduced the usefulness of 
data near 6 c . m , = 90°. However, a large enough range 
in 9 could be covered to demonstrate the change in the 
angular distributions with energy and the need for partial 
waves higher than Ss and Sp. The dashed curves show 
calculations taken from [pl| . 




















A, at 350 MeV 






A, at 375 MeV 



FIG. 4. The integrated pp — > pn7r + observables 
A<7T/cr un p i = -(A xx + A yy ), D = A xx - A yy , and A y as listed 
in Table 1 are shown as a function of rj, the maximum pion 
center-of-mass momentum in units of the pion mass. The 
observables A xx , A yy and A y were integrated over all kine- 
matic variables except the pion azimuthal angle. Solid curves 
are from Hanhart et al.. These calculations are for the full 
phase space and are not directly comparable with the data. 
The dashed curves approximate the theoretical calculations 
for the actual phase space observed. 

The integrated spin correlation coefficient S is related 
to the spin-dependent cross section by Aa' T /a to t = — S, 
and we find that these two analysis methods agree. The 
integrated analyzing power is given by A y = (— ae • 
ax) 1 / 2 . The resulting polarization observables are shown 
in Fig. 4 and listed in Table 1. These values, but not 




A, at 400 Me' 
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FIG. 5. The measured pj3 — » pn7r + analyzing powers 
Ay(# c . m .,T) compared with fits (solid curves) described in 
the text. The angular distributions suggest some asymme- 
try around 9 c .m. = 90°. Acceptable fits require the inclu- 
sion of an associated Legendre polynomial with = 2. 
The detector size did not permit reliable measurements for 
—0.25 < cos(0 c .m.) < 0.25. The dashed curves are calcula- 
tions by Hanhart et al.[18]. 
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The solid curves in Fig. 5 represent fits with the rela- 
tion 



(X)) * (a Pl(x) + b Pi(x) + c p£(x)), 

(3) 

where the expressions P 1 u (x) are associated Legendre 
polynomials and x = cos (6). Pj(x) was not needed or 
used for the fits shown; however, P^x), i.e. the inclusion 
of an L-K = 2 term, was essential to obtain a reasonable 
representation of the data at the four energies. 

We have found that the recent calculations of the 
Julich group are quite successful in predicting the inte- 
grated spin correlations for pp — > pmr + in spite of their 
difficulties with pp — > ppn . This dramatic difference in 
the otherwise similar reactions should help identify the 
reason for the less than satisfactory agreement of calcu- 
lations and measurements for pp — > ppir" . A further 
indication of remaining shortcomings is found in the dif- 
ference between calculated and measured angular distri- 
butions for A y (0 c . m .) in Fig. 5. It appears that some di- 
agrams that are secondary in pp — > pnn + but important 
in pp — > ppTr° have been omitted or treated inadequately. 

We are currently analyzing additional measurements of 
the pp — > pn7r + reaction with longitudinal beam and tar- 
get polarization. Such results together with the present 
data may allow a deduction of the dominant partial waves 
in a model independent way. 
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TABLE I. Product of beam and target polarization and the deduced integrated spin correlation coefficients. The observables 
listed are affected differently by the limited observed range for the polar angle 8. Model calculations show that the lack of 
counts for the range —0.4 < cos(# c .m.) < 0.4 makes the listed values for -S and D larger than what a complete integral would 
have produced. The effect is more pronounced for A y because A y (# c . m .) has a strong minimum near c . m . = 90°. The last 
column provides the full integral over the fits to A y (6) in Fig. 5 and a better value for A y . See text. 
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325.6 


0.464 


0.456±0.003 


0.553±0.038 


0.031±0.066 


-0.167±0.016 


-0.207± 0.019 


350.5 


0.623 


0.342±0.004 
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0.104±0.080 


-0.161±0.016 


-0.210± 0.020 


375.0 
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0.067±0.021 


-0.140±0.005 


-0.194± 0.016 


400.0 


0.871 


0.526±0.006 


1.274±0.019 


0.095±0.035 


-0.164±0.009 


-0.212±0.015 
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